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The asymptotics of some spin foam amplitudes for a quantum 4-simplex is known to display rapid 
oscillations whose frequency is the Regge action. In this note, we reformulate this result through 
a difference equation, asymptotically satisfied by these models, and whose semi-classical solutions 
are precisely the sine and the cosine of the Regge action. This equation is then interpreted as 
coming from the canonical quantization of a simple constraint in Regge calculus. This suggests 
to lift and generalize this constraint to the phase space of loop quantum gravity parametrized 
by twisted geometries. The result is a reformulation of the flat model for topological BF theory 
from the Hamiltonian perspective. The Wheeler-de-Witt equation in the spin network basis gives 
(3JT)' difference equations which are exactly recursion relations on the 15j-symbol. Moreover, the semi- 

classical limit is investigated using coherent states, and produces the expected results. It mimics the 
1 classical constraint with quantized areas, and for Regge geometries it reduces to the semi-classical 

| equation which has been introduced in the beginning. 
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INTRODUCTION 

Asymptotics of spin foam amplitudes from semi-classical Hamiltonian dynamics 



A good spin foam model for quantum gravity is (often) expected to reproduce Regge calculus (a large distance 
approximation of general relativity) in the classical limit. The idea goes back to Ponzano and Regge, [3], who 
\q , made the key observation that the Wigner 6j-symbol, an object from the theory of representations of SU(2), behaves 
7— I ' for large spins as the cosine of the Regge action for a tetrahedron, with the spins as edge lengths. This gave a 
, model for quantum gravity in three dimensions, where spins are interpreted as quantized lengths. The achievement 
of Loop quantum gravity (LQG) then gave a new birth and justification to the idea that quantum gravity can be 
■ formulated from algebraic objects, coming from the representation theory of a Lie group, attached to chunks of 
. . | spacetime (simplices, polyhedra). 

Still, it is not so obvious and straightforward to imagine why the semi-classical limit of spin foams would have to 
be expressed in terms of the Regge approximation to general relativity (remember spin foams are initially designed to 
provide transition amplitudes between the kinematical states of LQG, based on cylindrical functionals of the Ashtekar- 
Barbero connection, 17|). This idea was suggested in [19] . In particular, it is based on the fact that LQG supports 
a discrete area spectrum (built from the Casimir of SU(2)), very similar to the Ponzano-Regge ansatz for quantized 
lengths in three dimensions. 

So when a new model is proposed, the natural thing that is to be done is to check its semi-classical limit, where 
by "checking" it is usually meant chasing after the Regge action. However, it turns out that the models that have 
so far attracted the most attention all have such Regge contributions, in particular models that are known not to 
describe quantum gravity, like the Ooguri model (a model for the topological BF theory in four dimensions . 1 1 611 ) and 
the Barrett-Crane model (though the oscillatory part involving the Regge action is only a subleading term [13(). 

In a series of papers (see for instance 0, (gj), the asymptotics of the 15j-symbol (for the Ooguri model) and for 
the Euclidean and Lorentzian EPRL 4-simplex have been precisely studied. It appears that different behaviours are 
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observed according to whether or not the boundary data determine a Regge metric on the 4-simplex, so that oscillations 
with the Regge action do occur or do not. 

So we would like to get a criterion which would tell us whenever a model has such oscillations with the Regge action, 
say at the leading order. We obtain such a criterion as a difference equation of second order on the 4-simplex amplitude, 
([3]) . This equation is actually well-known from the three-dimensional case. Indeed the 6j-symbol is fully characterized 
by a second order recursion relation (coming from the Biendenharn-Elliott, or pentagon identity), which, although 



generally complicated, simplifies in the semi-classical limit where it allows to determine the asymptotics, [12l. \21\. 

Since the asymtotics provides a regime where the spin foam amplitude may be approximated by some quantum 
Regge calculus, it is natural to look for an interpretation, or even better, for a derivation of this difference equation as 
the quantization of a constraint in Regge calculus. It turned out to be very simple, and natural. The corresponding 
constraint is a sort of flatness constraint which enables to built a flat 4-simplex from its boundary. It states that the 
momenta conjugated to the triangle areas have to be the dihedral angles between adjacent tetrahedra, computed from 
the areas like in a flat 4-simplex. This gives a first link between a classical constraint and the asymptotics of spin 
foams. 

So far the relation between spin foam models and the Hamiltonian constraint of general relativity has been partic- 
ularly evasive (see open problem (14) in [20]). Our result gives a taster for this relation. To go further and ultimately 
savour it, we need to derive the semi-classical constraint from the quantization of a Hamiltonian operator in LQG. 
We will perform this task in the case of the Ooguri model. The analysis is a simple extension of results to appear 
from a collaboration with L. Freidel, There it is shown in the 3d case that a projection of the curvature onto the 
components of the triad, thus taking the form of the Hamiltonian constraint EEF, can be quantized in LQG. In the 
simplest situation, on the boundary of a tetrahedron, the Wheeler-de-Witt equation is a difference equation which is 
exactly the recursion relation defining the 6j-symbol. In 4d, the physical (flat) state on the boundary of a 4-simplex 
is the 15j-symbol. One can lift the Hamiltonian used in 3d to 4d, and using the methods of [§[, we claim that the 
Wheeler-de-Witt equation reproduces the recursion relations satisfied by the 15j-symbol which were derived in 

The organization is as follows. In the section [IJ 

• we exhibit a difference equation whose solutions are the exponentials of ±i times the Regge action of the 4- 
simplex. 

• we derive this equation as the quantization of a classical constraint in area Regge calculus. The constraint states 
that a point on the phase space is given by a set of ten areas and ten dihedral angles which are those of a flat 
4-simplex, determined by the areas. 

The section HT1 focuses on the Wheeler-de-Witt equation for the Ooguri model. 

• We define a classical constraint, attached to a node and a cycle of a spin network graph, by projecting the 
curvature onto some components of the gravitational field, III Al 



• This Hamiltonian is rewritten in terms of twisted geometries [14j . a nice parametrization of the LQG phase space 
on a single graph. It appears as a generalization of the above constraint for Regge calculus to the whole LQG 
phase space (including non- Regge boundary data). |HB] 

• The corresponding Wheeler-de-Witt equation is studied, more particularly in the large spin limit, in the coherent 
state basis. Using the WKB approximation, it reduces to the classical Hamiltonian on twisted geometries, with 
quantized areas. 

In particular, in the Regge sector of boundary data, it reproduces the semi-classical equation asymptotically satisfied by 
spin foam models. The natural variables are areas and normals of triangles. Thus, it strengthens from the Hamiltonian 
point of view the result that quantum area-angle calculus is the semi-classical limit of quantized geometries (in the 
Regge sector). 
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We will also argue that such difference equations obtained through canonical quantization in the LQG framework 
lead to a position where the same analysis as that of Q can be done and used to extract the asymptotics (like in 3d 



a pc 

actually, @,l2l|) 



In the section UTil we discuss additional interesting difference equations. One is derived from our main semi-classical 
equation ([3]) , and shown to probe the closure of the simplex (it was already introduced and precisely described in Q , 
though from a quite different path). We also sketch the possibility of introducing more speculative constraints, whose 
asymptotical behaviour exhibits oscillations with the Regge action. 

All technical details are skipped in the main text to ease a fluent reading, and are reported in appendix. 



I. A NEW LOOK AT THE SEMI-CLASSICAL BEHAVIOUR OF 4-SIMPLEX SPIN FOAM 

AMPLITUDES 

Crucial references on the asymptotics of spin foam models are [3- 6] . There it is shown that several spin foam models 
get in the large area limit rapid oscillations with the Regge action. We want to track back this phenomenon to the 
fact that they satisfy in this regime the same equation, solved by exponentials of i times the Regge action. 

Consider a 4-simplex, with five tetrahedra on its boundary labelled by a = 1,...,5. The semi-classical regime 
corresponds to large values of the quantum numbers of triangle areas (j a b S y)a<t>- A spin foam amplitude for the 
4-simplex is determined by some boundary data, including these quantum areas. The additional boundary data are 
coherent states labelled by spinors in the Lorentzian EPRL model 0], and by points on the 2-sphere (N a b G S 2 ) a ,b 
for the SU(2) Ooguri model Q|. The latter have been introduced in spin foam models to gain geometric control on the 
quantum amplitude, . Indeed, there are some distinguished sets of boundary data which allow to construct the 
geometry of genuine flat 4-simplices in 1R 4 . There, the vectors (N a b) are the normals to the triangles with areas j a b, 
and the five tetrahedra consistently glue to form a 4-simplex. This is the geometric sector of the model, consisting 
in all possible Regge geometries. In this sector, it has been shown in a series of paper (see Q and other references 
therein) that for homogeneously scaled areas (Xj a b), with A ^> 1, 

W(\j ab ) = N+ e lAS *(^) + e - ixSR< - j °- b) + o(l), (1) 

up to a global scaling. The coefficients N± depend on the boundary data, but not on A. The Regge action is as usual: 

a<b 

where the areas are actually given by: A a b = j a b in the Ooguri model and A a b — j j a b in the EPRL model, the global 
parameter 7 being known as the Immirzi parameter (this last point is in agreement with the area spectrum of canonical 
Loop Quantum Gravity) . The angles Q a b are the dihedral angles determined by the geometry of the boundary data, 
which arc then the angles between the boundary tetrahedra. 

So the large spin regime provides a notion of semi-classical quantum gravity, where the spin foam amplitude is 
well approximated by some quantum Regge calculus. Still, because this happens for several models in spite of their 
fundamental differences, one would like to understand this behaviour in a somewhat universal way. This is what we do 
by exhibiting a difference equation asymptotically satisfied by these models, which directly leads to this semi-classical 
approximation. Then we derive this equation as a naive quantization of some flatness constraint in area Regge calculus. 
To go further along these lines, it is desirable to get an equation satisfied by the full amplitude, and not only in the 
semi-classical limit. 



1 There are also some phase ambiguities, which will play an important role in the next sections, but we do not need such details at this 
stage. 
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A. Recursion relations in the asymptotics 



Our first point is that linear combinations of exponentials of ±i times the Regge action are the solutions to the 
following difference equation of the second order: 



A ab + 2 (1 - cos6 a6 ) V(A ab ) = 



(3) 



when solved via the semi-classical approximation for large A. Here A ab is the discrete second derivative with respect 
to the area variable A ab (may it be j ab or jj ab ): A/(x) = f(x + 1) + f[x — 1) — 2f(x). An equivalent form (which will 
be that naturally coming out in the next sections) is obtained by defining some ladder operators which shift an area 
by±l, 



$i b V(A ab ) = V(A ab + 1), S- b V(A ab ) = V(A ab - 1). 
Then, the semi-classical equation becomes: 



(4) 



2( 6 tb+ S ab) - COS ©ah 



V(A ab ) = 0. 



(5) 



We look for solving the equation a la WKB, when all spins are rescaled by A ^> 1 , and with the ansatz: 

^(\A ab )=<f>(A ab )e iS ^ A ">\ (6) 

We assume $ does not scale with A, while S scales linearly, so the idea is as usual: a slowly varying amplitude, with 
a rapidly oscillatory phase. To zeroth order, 

i>(\A ab ± 1) ~ ip(\A ab ) e ±i3 '^\ (7) 

where S' is the derivative of S seen as a function on the real line. The equation ^ becomes: 

cos S'(A ab ) - cos 6 a6 = 0, (8) 

or: S'(A ab ) = ±Q ab . So one has to integrate the dihedral angle with respect to the area. The result is known to be 
the Regge action, since when varying it with respect to A ab , the variations of the dihedral angles cancel thanks to the 
Schlaefh identity, Y, a <b A <A S ®ab = 00. 

To better control the approximation made in (J3]), it is useful to derive it directly from the asymptotic analysis of 
the models. The vertex amplitude W is exactly defined as: 

W= [dfi(X) exp(-^ J 4 cd s cd (A)). (9) 

X denotes the set of all variables to be integrated, and the function s ab has a positive real part. Thus, it can be 
evaluated using a saddle point and stationary phase method when rescaling the spins by a large A. Applying the 
operator A a fc gives without approximation: 

A ab W = 2 { d M (l)e-S^W (cos(s afc (A)) -f). (fO) 



2 Strictly speaking, the Schlaefli identity is usually considered for any variations of the edge lengths. In the geometric sector, the set of 
spins and coherent states enables to reconstruct these lengths. Thus, the chain rule allows to write the Schlaefli identity for variations 
of the areas also. 
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In the large spin limit, we can just evaluate cos(s a &(X)) on each saddle point. If it is the same for all of them, then 
it can be factorized from the amplitude. This is what happens in the SU(2) Ooguri model and the Lorentzian EPRL 
model, where the saddle points X* gives: s a b(X*~) = ±O a t>. This point is important and prevents the asymptotics from 
getting other frequencies than the Regge action itself. As an example, our reasoning does not apply to the geometric 
sector of the Euclidean EPRL model, since it receives oscillations from the Regge action and from 7 -1 Sr also, with 
the same scaling. Then, the spin foam amplitude satisfies a higher order difference equation in the asymptotics, which 
is simply the product of the difference operator Q for both frequencies: 

V(A ab ) = 0. (11) 



A„ 



2 1 



COS ) 

'V 



2 (l - cos9 q6 



B. Recursion relations as Wheeler-de-Witt equations in quantum area Regge calculus 

Since semi-classical spin foams can be approximated with quantum Regge calculus, we now want to understand our 
main equation ^ in this framework. Here comes our second important point: this equation Q has a nice geometric 
interpretation as a quantization of a constraint in Regge calculus. Consider the set of areas (A a i,) such that they 
uniquely determine a genuine flat 4-simplex as the configuration space. Like in [Io| . we take the conjugated momenta 
to be angles {Q a b)i with the canonical brackets: 

{A a b, Ocd} — S{ab),(cd)- (12) 

This framework has been derived from a canonical discretization of the Plebanski's action for gravity in [ljj (see 
there for the full details on the symplectic structure of the phase space). On this phase space, we also consider the 
constraints: 

Xab = COsdab - COsQab(A) = 0. (13) 

This was already studied in the above reference, and there argued to form an Abelian algebra. It should be noted 
that the authors of fici| were then interested in the gauge symmetry corresponding to the translation of a vertex of 
the simplex. Here, we would like instead to generate independent shifts of areas to produce our equation of interest 
([3]). For that purpose, the constraint (j!3[) is what we need. Furthermore, its geometric meaning is quite clear: the 
momenta (9 a b) are constrained to be the dihedral angles (Q a b) of the flat 4-simplex determined by its areas. 

Let us now proceed to the most naive quantization, using wave functions of the angles. They can be expanded onto 
the Fourier components, (e* £-?<»<> 9 <j'') ) where the integers (j a b) are the eigenvalues of the area operators A a b (they get 
discrete spectra since the variables 6 a b live on a compact set). Then, periodic functions over (0 a b) act by multiplication, 
and in particular: 

[Y, Hjcd)e 1 ^^] = ]T Mjab T UcdK^'H (14) 

{jed} {jed} 

= 4>- (is) 

This simply means that on the Fourier coefficients ip(j c d) of a state \ip), the operator e ±l9ab acts by shifting the 
variable jab by =Fl. Also, as the Fourier exponentials are the eigenfunctions of the area operators, we simply promote 
the complicated functions Q a b(A) to operators through: 

cose^A)^^* = cos 9 ab {j) i " d9 - d , (16) 

as far as the set {jed) allows to define the dihedral angles. Thus the classical constraint \ab can be imposed at the 
quantum level, 

X^)=0, (17) 
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where it becomes exactly the difference equation we are looking for: 

A afc + 2 (1 - cos e a6 )] Tjj(jab) = 0. (18) 

II. QUANTUM DYNAMICS OF THE FLAT 4-SIMPLEX IN LOOP QUANTUM GRAVITY 

(THE OOGURI MODEL REVISITED) 

The above approach has obvious limitations: 

• First, it only holds asymptotically, and we expect the full quantum gravity amplitude to satisfy a difference 
equation with non-trivial coefficients, which would contain all information about the full asymptotic expansion. 

• Second, it is not clear what the role of the additional boundary data of spin foams (the normals to the triangles) 
can be here. However, they are part of the phase space of Loop Quantum Gravity (on a single graph). In 
addition, it has been shown in previous studies that spin foams are better understood in terms of area-angle 
Regge calculus, 0], instead of area calculus (well-known to suffer from several drawbacks), 11 1. This leads to 
the third point. 

• We have so far focused only on Regge geometries, since the asymptotic behaviour is different on the other 
configurations. But from the LQG point of view, there is no specific reason to distinguish between Regge 
and non-Regge geometries. Furthermore, the Ooguri model, which also shows up these different asymptotic 
behaviours, is nevertheless built from a single constraint, namely the flatness of a gauge field, without regards 
for the amount of geometricity contained in the canonical momenta. 

This leads us to revisit the Ooguri model for SU(2) BF theory, with a new form, more geometric, of the Hamiltonian 
constraint. 



A. The proposal: projecting the curvature 

So we now turn to the phase space inherited from LQG on the dual complex T to the boundary of a 4-simplex. 
The notation a = 1, . . . , 5 is kept for tetrahedra of the triangulation, and hence (ab) for triangles and (abc) for edges. 
By duality, they correspond on the complex T to nodes a = 1, . . . , 5, links (ab), and faces, also referred to as cycles 
(abc). We will mainly use the terminology corresponding to cells of T, but also switch to the point of view of the 
triangulation as soon as we find it relevant. 

The phase space is precisely the same as that of SU(2) Yang-Mills on such discretization, and the same as in the 
topological SU(2) BF theory, up to a scaling of the fundamental brackets by the Immirzi parameter (that we will ignore 
in the following). The ten links of T carry SU(2) elements, (g a b), which represent parallel transport operators between 
the source and target vertices of each link|f|. The phase space is then simply the cotangent bundle over SU(2) L , with 
L = 10 and with its natural symplectic structure. The momentum to g a b is thus a 3-vector E a b, for a <b, with 

{Kb,9ab} = T*g ab . (19) 

These momenta can be seen as smearings the triad field Ef of the continuum over the triangles (dual to the link of T) , 
and will be therefore called triad variables. We have denoted (r*)i = i : 2,3 anti-hermitian generators of the Lie algebra. 
The standard interpretation is the following. Each tetrahedra of the boundary carries a local reference frame, and E a b 



3 We choose the notation so that g a f, goes from b to a. Also: g^ a = g £ . 
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is defined relatively to that of the tetrahedron a. Momenta E ba acting on the right of g ab can be defined by parallelly 
transporting E ab to the frame of the tetrahedron b using the adjoint action of the group: 

E ba = -Ad(gJ)E ab . (20) 

Like in lattice gauge theory, the gauge group is SU(2) V , where V = 5 is the number of nodes of the graph. 
We now consider a Hamiltonian constraint for the Ooguri model. It is usually taken to be: 

9(abc) = 9ab9bc9ca = 1, (21) 

That is to say: the parallel transport around the cycle (abc) made of the three links (ab), (be), (ca) is trivial. But it 
is unlikely that its quantization will help to understand quantum gravity. So we would like a constraint which would 
look like closer to the Hamiltonian constraint of general relativity, e 3 k EfE^F(A) k a ^ = (where a, (3 are indices on 
the canonical surface, and A is a SU(2) gauge field of curvature F and conjugated momentum E). First the curvature 
is discretized around the 2d regions of the complex dual to the triangulation. More precisely, consider the region 
bounded by the links (ab), (be), (ca), then the component of F along these directions is regularized in LQG as: 

513 - (Ad(<7(a6c))f ■ (22) 

Then, the idea is to project it at each node of the cycle along the two triad variables which meet there. Define: 

Hbc — ^ab ■ E ac — E ab ■ KA(g^ abc ))E ac , (23) 

and the constraint H£ c — 0. For each cycle, there are three such constraints (and for a generic triangulation, there 
is one constraint for each node of a cycle). So there are enough constraints to enforce <?( a f> c ) = ±1 when the triad 
variables around the cycle span the three dimensions 0- 



B. A Hamiltonian for twisted geometries 

We are interested in the relation between H bc and the previous constraint \ab, (fT3|) . at the classical level, and 
with the semi-classical equation ([3]) after quantization. Classically, it is possible to define in the geometric sector the 
dihedral angles Q a b via the triad variables, and like in (icj |. another notion of dihedral angles involving the group 
elements. Then, the above constraint just states the equality of the two notions. Details will appear in a collaboration 
with L. Freidel 0. 

Here, we prefer to translate the constraint into the language of twisted geometry [14(. This is a nice reparametrization 
of the LQG phase space which makes clear the nature of the involved geometries. In particular, space is formed by 
genuine polyhedra like tetrahedra, but their gluing does not lead to Regge metrics, since two adjacent polyhedra may 
describe their common boundary with different shapes. We take advantage of this fact and give an interpretation of 
H bc which also holds for non-Regge situations. 

The parametrization maps the set (E ab ,g ab ) to a new set: 

(E ab , gab) -> (A ab ,N ab ,N ba ,£ab), (24) 

defined as follows. A ab is the norm of E ab (and equals that of E ba ), and N ab its direction: 

E ab = A ab N ab , (25) 



4 At least, it is not hard to see that there is no smooth deformation of this relation. However, there may be a finite set of possibilities that we 
have not investigated. In particular, the solutions of the equation E\ -E^ — E\ ■ Ad(g)i?2 = are: g = exp(ii_Ei) exp(r;(_Ei X E2)) exp^^): 
where ti, t2 are arbitrary. But r] admits only a finite number of values, since there is a finite number of SO(3) rotations with axis E\ X E2 
solving the equation. 
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for all a, b. Then, since N a b and Nba are taken as independent, the equation (|2U)) : N a b = — Ad(g a b)Nba has to be 
solved for g a b- Take a set of SU(2) rotations (n a b(N)) a .b such that n a b maps an axis of reference in R , say z, onto 
N a b. This leads to the introduction of the angles £ a f, through: 

9ab = n ab e e^ abTz n^. (26) 

(The matrix e= (-10) ^ there to account for the minus sign in the parallel transport relation (|20j) . since e maps the 
direction z onto its opposite —z.) Obviously, the normals and the triad variables are unchanged when adding a phase 
on the right of n a b like: 

n ab -> n ab e x ^. (27) 

The invariance of g a b then requires to change t; a b accordingly. In particular, for a given g a b, £ a b can always be 
reabsorbed into the rotation n a b or rib a - The set of rotations (n a b) is very convenient, as we will see, and its use 
prefigures what happens at the quantum level. Indeed, the semi-classical coherent states we will later use are actually 
labelled by such rotations rather than only by the normals (N a b)- 

The generator of gauge transformations on the vertex dual to the tetrahedron a is: 

Aab N ab = 0, (28) 

This condition actually takes the form of a closure relation for the tetrahedron and hence leads to this nice interpre- 
tation: the variable A a b is the area of the triangle (ab), while N a b and Nba are respectively the normals to the same 
triangle with respect to the frame of the tetrahedra a and b. So it guarantees that one can built a flat tetrahedron in 
R 3 for each a. 

The areas and normals describe the intrinsic geometry of the canonical surface. In particular, in the gauge invariant 
sector, the dihedral angle (f>^ c between the triangles (ab), (ac) is given by: 

1 a E a b ' E ac /on\ 

cos(p bc = — - — — . (29) 

siab s±ac 

Another key quantity we will need, defined only in terms of the six normals around a cycle (abc), is: 

bc K ' sin^sin^ V ' 

If this quantity is independent of a (that is computing it from any cycle containing the link (bc) gives the same answer) , 
then it is exactly the 4d angle Qbc between the tetrahedra b and c, computed from the normals. This is exactly the 
criterion that turns a set of variables satisfying (|28p into a Regge metric on the triangulation [ll| : 

cos 6^ (A?) =cose^' ) (iV). (31) 
The splitting between intrinsic and extrinsic geometries in the twisted parametrization, and in particular, the 



information about the extrinsic geometry contained in the set of normals (N a b) has been discussed in Here, 

we are able to go further on this issue by writing the constraint H£ c in this new set of variables. To get a definite 
expression, we need the Euler decomposition of the product n~£n ac : 

Kb^ac = e a ^ T ' e<*-*?>« e <* T » , (32) 

which defines the angles (a bc ) a ,b,c (notice however that they are changed under ([27])). 

The Hamiltonian H bc then admits the following form on twisted geometries (see appendix): 

Hb c = -A ab Aac (COSC - COS 4> a c COS <j>ba + sin 4>ac ^4>ba C0S (6c + «L + "L)V ( 33 ) 
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The surprise is that this is the form of the standard relation between the 3d and 4d dihedral angles within a flat 
4-simplcx, (f3T))) . though it holds on the whole phase space. Now restrict attention to the Regge-geometric sector 
(where (f3"Tj) is satisfied). As soon as the tetrahedra are non-degenerate (the 3d angles are neither nor ir), the angle 
(6hc + a ta + a ba) can be extracted from the constraint H£ c = 0, to give the dihedral angle <d bc . Moreover, since £t, c 
is also independent of a obviously, then {a b ca + a£ Q ) also is. So we can use (|2"Tf to achieve a phase choice where this 
combination is ir. This gives 

H£ c oc cos Q bc (N) - cos(£ hc ) = 0. (34) 

So the classical constraint H£ c = really corresponds to the constraint \ ab = in their common domain of applicability 
(the geometric sector), and is clearly related to building a flat 4-simplex out of its boundary tetrahedra. A key difference 
is that the dihedral angles of the 4-simplex are rather computed from the set of normals {N ab ) rather than from the 
areas. This suggests that the semi-classical limit of LQG is given by quantum area-angle Regge calculus. Further, it 
makes the formula applicable on the whole phase space, by going back to (133|) for H§ c , since the latter still makes sense 
outside of the geometric sector (when the gluing of the tetrahedra is not that of a 4-simplex), and also for degenerate 
tetrahedra (when some 3d angles are such that sin0£ c = 0). 

We would like to mention that at this stage it is possible to discuss the solutions of the constraints (|3"3"]l in a way 
which is fully parallel to the analysis of [§]. First, in the geometric sector, there are clearly two solutions, 

6 C = ±e bc (A). (35) 

Then, assume there are at least two distinct solutions, (£ bc ,ti bc )- That leads to: E ba ■ Ad(<7^,)£' ca = E ba ■ Ad(g^ c )E ca . 
This relation has been studied in [loj where it was coined edge-simplicity constraint, and shown to actually imply the 
Regge gluing relations (|3"Tj). So in the non-geometric sector, the constraint has either one solution, or no solution. 



C. The Wheeler-de-Witt equation and its semi-classical regime 

1. The Wheeler- de-Witt equation as recursion relations 

When a theory has gauge symmetries, the latter turn into constraints in the Hamiltonian analysis. They can 
be imposed either before or after quantizing. In general relativity and BF theory, the Hamiltonian itself is a 
constraint, and the program of LQG is to quantize first and then constrain. The kinematical Hilbert space is 
Uv = L 2 (SU(2) i /SU(2) y ), spanned by the so-called spin network functions. These are just built from the Fourier 
modes of the ten group elements, that are their matrix elements in the representation (j ab )a< b , while the magnetic 
numbers are contracted with a specific tensor (t a ) on each node a, called intertwiner, which ensures gauge invariance. 

s W4 (5ab)= Y[{j ab m ab \g ab \j ab m ba )l[ii m ^. (36) 

{rn ab } a<b a 

The quantization of operators is direct. Gauge invariant functions of the group elements act by multiplication, while 
E % ah acts as a left derivative: 

E lb (jab m ab \g ab \j ab m ba ) = i(j ab m ab \T l g ab \j ab m ba ). (37) 

In particular, the square of the above equation produces the area spectrum, that of the Casimir of SU(2). 

The basic building block attached to every 4-simplex in Ooguri's model for SU(2) BF theory is a Wigner 15j-symbol. 
The latter can also be seen in the Hamiltonian framework as the flat state on the boundary of a 4-simplex, satisfying 
(|2Tj) on each cycle, expanded in the spin network basis. 

Though it is not hard to express the action of the operator E ab ■ E ac on a spin network state, the result actually 
depends on the choice of a basis of intertwiner. A standard basis is obtained at each node of the graph by pairing 
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the links meeting there, expanding the tensor products of their representations into irreducible representations, and 
choosing a spin common to the two tensor products. If (ab) and (ac) are paired together, then the operator E ab ■ E ac 
is diagonal on spin network states, Q. As for the operator E ab ■ Ad(<7( a ;, c ))i? ac , one first rewrite it as: 



E ab ■ Ad(g( abc ))E ac = E ba ■ Ad(g bc )E c 



(38) 



The triad variables produce insertions of generators on the link (ab) at the node b, and on the link (ac) at the node 
c. The action is not diagonal since Ad(<?& c ) is a multiplication by matrix elements in the spin 1 representation. This 
produces shifts on the spin j bc to j bc + k, for k = —1, 0, 1. 

Depending on the chosen pairings on the spin network nodes, the action of H c b l c thus leads to different equations (all 
of them being difference equations of the second order acting on one or more spins). However, using the same tools and 
methods as those of [8j , it can be oberved that these difference equations exactly take the form of recursion relations 
satisfied by the Wigner SU(2) 15j-symbol @. On the one hand, this 15j-symbol is the basic building block attached 
to every 4-simplex in Ooguri's model for SU(2) BF theory. On the other hand, it is in the Hamiltonian framework 
the flat state on the boundary of a 4-simplex, satisfying (I2ip on each cycle, when expanded in the spin network basis. 
The recursion relations have been derived in from the topological invariance spin foam model - a method which 
already suggested a strong connection to the classical symmetry of the theory. 

Assume for instance that the intertwiner i\ pairs (712,713) together to a virtual spin i 1; that l 2 pairs (712,723) to i 2 
and 63 pairs (723,713) to a virtual spin 13. Then, the Wheeler-de-Witt equation: 



#23 \f) = 0, 

is actually the same as the recursion relation on the 6j-symbol, 

A-1U12) 4>(ji2 - 1) + AoUiz) ip(ji 2 ) + A +1 (j 12 ) ip(j 12 + 1) = 0. 
The coefficients A±i take the form: A + i(j) = jE(j + 1), and A-i(j) = (7 — l)E(j), for 

E( Jl2 ) = [((il3 + h + I) 2 - Jl 2 2 ) (fl ~ (J13 - *l) 2 ) ((J23 + * 2 + I) 2 ~ Jl 2 2 ) (Jl 2 2 - 0*23 - *2) 2 ) 

and the coefficient Aq is given by: 



(39) 



(40) 



(41) 



AoO'12) = (2ji2 + 1){2[j 13 (j 13 + l)i 2 (i 2 + 1) + 733(733 + l)*i(*i + 1) - 712(712 + l)ia(*3 + 1)] 

- [713(713 + 1) + »i(*i + 1) - 712(712 + 1)] [723(723 + 1) + i 2 (t 2 + 1) - 712(712 + 1)] }• (42) 

If now on the node a = 1, the intertwiner i a pairs (712,715) together to the virtual spin i\, then the equation becomes 
more complicated: 



^2 A ei2tei (j 12 ,h)ip(ji2 + £12, ii + £1) = 0. 



(43) 



ei2,ei = -l,0,l 



Such relations were derived by writing down explicitly a special invariance of the Ooguri model under a change of 
triangulation. Thus, it was known that these relations encode the symmetries of the model at the quantum level. But 
they had so far never been derived from the quantization of a Hamiltonian constraint. 



The generic process is that the action of a triad variable inserts a generator, and then, the contraction of their vector indices produce 
graspings on the spin network in the spin 1 representation. After some recoupling, one can extract a special 6j-symbol with a spin 1 at 
each node of the cycle (including the virtual spins of intertwiners). 
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2. The Wheeler- de- Witt equation in the semi-classical regime 



The above equation is not really suitable for the semi-classical analysis however. It comes from the uncertainty 
principle, that a tetrahedron is described quantum mechanically by only five quantum numbers (four areas (jab), and 
one spin i a to specify the intertwiner) 0. So to launch our Wheeler-de-Witt equation in the semi-classical limit, we 
go to an overcomplete basis of coherent intertwiners [l5j]. First build the usual SU(2) coherent state \j,n(N)) from a 
SU(2) rotation n(N) which maps the reference axis z onto a unit 3- vector N: 

\j,n(N))=n(N)\j,j). (44) 

It is important to keep in mind that the state is not fully determined by the direction N, but also by a choice of 
phase. Indeed, changing n like in (|27p does not affect the vector TV, but multiplies the state by a phase. A coherent 
intertwiner L a (n a b) on the tetrahedron a is labelled by four rotations (n a b)b^a corresponding to four unit vectors of 
K 3 , {N a b)b^a- It is defined by a group averaging process: 

dh a (g>b=£ah a \j a b,n ab ). (45) 

SU(2) 

It is shown in [15] that the norm of this intertwiner is peaked for large spins on vectors (N a b) which satisfy the closure 
condition (|28[) . but with quantum areas (j ab )- Therefore, these vectors can be interpreted as normals to the triangles 
of a (up to a global rotation). The spin network state in the basis of coherent intertwiners reads: 



s {j a b,n ab ] 



5 

(9ab) = / TTc#ia \\{3ab,n ab e\h a 1 g a bh b \j ab ,nba)- (46) 



(The matrix n a bt sends the axis of reference z onto —N a b- This choice is a matter of convenience.) 

To evaluate the action of Hjf on such states, some work is necessary, which is reported in appendix. However, the 
result is very natural. Indeed, the data of these coherent states match the data (A a b, N a b, Nb a ) of twisted geometries. 
So it is expected, if they really have a nice semi-classical behaviour, that they will lead to some simple quantum version 
of the classical Hamiltonian on twisted geometries (|33|) . 

The nice semi-classical behaviour is inherited from the following property of SU(2) coherent states: 

r\j,n(N)) = (-i)jN\j,n(N))+o(j). (47) 
Since the triad operator E a b inserts a generator like this, we can guess the action of E a b ■ E ac : 

^E^TEac) fltt-»'»-»> ~ JabJac COS C fl ««*-«-»}. (48) 

Acting with Eb a ■ Ad(gh c )i? ca is a bit more involved, but there are no conceptual difficulties. The adjoint action of 
gbc is a Wigner matrix in the representation of spin 1, which recouples with the matrix elements of §b c in the state. 
This produces some shifts of the spin jbc to jbc + f], for 77 = —1,0,1. These shifts extend to the coherent states, to 
\jbc + fifcc), (jbc + ?7, n c bc\ in the large spin limit, via the operators Eb a , E ca . We thus get: 

(E ba ■ Ad(gbc)E ca ) S ^ n ^ ~ cosC cos^ a - i S in0 ac sin^ a (V l(<+aL) S+ + e i ^+ a ^ 8^ >=, (49) 
where the angles a are determined by: 

n ab n ac = e <cTl e (7r -^ )r » e<» T * , (50) 



like for classical twisted geometries. 
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The final step is to rewrite the quantum condition: 

Hg c W = 0, (51) 

for the coefficients of an arbitrary expansion in the coherent spin network basis. Again, we have to invoke the semi- 
classical, large spin limit, since coherent intertwiners generally have non-trivial overlap. They become orthogonal for 
large spins, and then the Wheeler-de-Witt equation for H bc reads: 

(cosC - cosC cos^ a ) i,(j bc ) + X - sin^ c sin <p c ba ( e *C««+»?. tyUc + 1) + e-^+^VO'&c - 1)) = 0. (52) 

Here the dependence on other variables than j bc have been dropped. This is our key equation, which generalizes 
the semi-classical equation satisfied by the exponential of the Regge action Q to the whole phase space of twisted 
geometries (with sufficiently large spins). 

The following is devoted to specializing (|52t to Regge or non-Regge boundary data, and see that it reproduces the 
equations and results previously discussed. We look for solving the equation a la WKB, when all spins are rescaled 
by A ^> 1, and with the ansatz: 

^(Aj fec ) = $(i6cK S(Aj H (53) 

We assume $ does not scale with A, while S grows linearly. To zeroth order, ip(Xj bc ± 1) ~ ip(\j bc ) e ±tS ( A:,6c ), where 
S' is the derivative of S seen as a function on the real line. So, the equation (j52"j) becomes: 

(cos (j% c - cos 4> b ac cos 0g a ) + sin 4> b ac sin <jf ba cos(S'(j bc ) + a b ca + a c ba ) = 0. (54) 

We can now make contact with the first semi-classical equation of the paper via the classical constraint H bc written 
for twisted geometries in (13"3"|) . Indeed, the latter is the same as this, with S'{j bc ) instead of the angle £ bc . 

Assume non-degeneracy of the tetrahedra, so that the formula (f30|) for the 4d dihedral angles <d ab as functions 
of the normals (N ab ) is well-defined. Further assume the boundary data satisfy the gluing constraints ([3"Tj) and the 
closure relation (|28|) . Then, we know from the discussion above (|34|) that: (i) the quantity into brackets on the left of 
([54]) is: sin0^ c sin<f) ba cos6;, c , (ii) the angle {a b ca + a ba ) can be set to tt by a change of phase in the coherent states 
\jab,n ab (N)). Then, we get to ©, 

cos S'(j bc ) - cos 8 6c = 0, (55) 

solved by exponentials of ±i times the Regge action S*r. Now assume that (|54|) has only one solution for S'. The 
same way as it is discussed in [3j , the value of S' can be reabsorbed into a change of phase of the coherent states (this 
change of phase depends on the boundary data {jab,N ab ,N ba )). Such a choice cancels the oscillations on the wave 
function, S = 0, like in the asymptotic analysis of [6j. 

III. OUTLOOK 

The equation ([3]) has been shown to come from a constraint stating that the momenta conjugated to the areas have 
to be the dihedral angles of a flat 4-simplex with these values of area (in the geometric sector). So this equation 
encodes the full information on the reconstructed 4-simplex, and we can imagine deriving from it other equations 
containing a few less geometric content, but still relevant. We here give one interesting example, which is a difference 
equation probing the closure of the simplex. Indeed, we can form the 5x5 Gram matrix: 



(Gab) = (cosG o6 ). 



(56) 
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Since the angles are those of a flat 4-simplex (actually determined by the areas and normals (A a },, N a b, Nb a )a<b, but 
this information will be lost in the coming equation), its determinant vanishes: 

det (G ab ) = 0. (57) 

From ([3]), we know that semi-classically, the multiplication by cosO a b can be compensated with shifts on the spin, 
by ($ab + ^ab)- This leads to the definition of a closure operator which annihilates the 4-simplex spin foam amplitude 
V(jab) in the large spin limit: 

det (\(S++S- b )) V(jab)^0. (58) 

This equation also holds for the Lorentzian EPRL model with shifts of the areas -fj a b instead. It can be directly 
checked, and this has been done in that the exponential of i times the Regge action is indeed a solution in the 
large spin limit, the key technical point being the use of (l57l) . 

It is remarkable that a 4-simplex spin foam amplitude satisfying this recursion relation in an exact way is actually 
known: the Barrett-Crane lOj-symbol (see [9| jl. Now that we understand that (|58l) is naturally a consequence of the 
flatness constraint (or semi-classical difference equation dU)), this result for the lOj-symbol is even more interesting. 
Indeed, it is unlikely that the lOj-symbol is annihilated by the operator [^(5 + + — cos0], since the rapid oscillations 
due to the Regge action only appear in subleading orders in the asymptotics. Nevertheless, it turns out that it satisfies 
the weaker constraint corresponding to the closure. 

The Hamiltonian H^ c we defined for the Ooguri model can be seen as a discretization of a quantity, e 3 \Ef -F 1 ^ 
for a fixed a and /? (that is a fixed face of the complex dual to the triangulation). So it is natural to try to implement 
the sum of the space indices to realize the Hamiltonian constraint of general relativity. A simple (naive ?) idea is 
to sum over the faces (cycles) that meet at the vertex a on the boundary of a 3-cell, say (abed). (Viewed from the 
triangulation, it means summing over the edges which meet a node within a tetrahedron), 

HI + H a od + H a db = 0. (59) 

From our analysis, it is clear that exponentials of ±i times the Regge action are annihilated in the large spin regime by 
this constraint (the three contributions would cancel independently). But we expect that there may be more solutions 
than for H^ c = which leads to the 15j-symbol at the quantum level. Since generates in some way the symmetry 
and the moves at the core of the topological invariance of the Ooguri model, the above constraint would relax that to 
asking for an invariance under a combination of the moves. This idea was argued in [9] on a toy example, but we may 
consider the above quantity as a more interesting realization. 

Such an operator is labelled by a node and 3-cell (instead of a node and a cycle for Ht c ). Note that in the case 
we have considered, that of the boundary of a 4-simplex, these 3-cells are actually of tetrahedral form. So pictorially 
speaking, this strongly echoes the recent proposal [l| of regularization of the Hamiltonian constraint in LQG. The 
latter consists in adding a small tetrahedral graph at the node of a spin network grah. An important feature is that 
a 15j-symbol could be extracted there. So it is definitely a good direction to compare the simple approach presented 
here (which perfectly works for the topological Ooguri model) with the standard approach based on the Thiemann 
trick, and more specifically on the volume operator of LQG (the latter does not seem to play any role here, a priori). 
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Note that the natural areas in the Barrett-Crane model are (2j a t, + 1), and not the spins themselves. So the operators 8 ab should then 
be understood as generating shifts of ±i in the spins (or equivalently, shifts of ±1 on the areas). 
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APPENDIX: TECHNICAL DETAILS 
On the Hamiltonian for twisted geometries 

Since the normals N ab can be defined by: 

N ab = Ad(n ab ) z, (60) 
we can easily use the rotations (n a b) to capture the geometric information. Let us form: 

n- h x n ac = e a ^ e (»-«J^ e «>». ( 6 1) 

This is the Euler decomposition of the product n~ b n ac . All that we need can be extracted from the matrix elements 
of this product in the representation of spin 1. Indeed, one can check that 0£ c as it appears above is the (3d) dihedral 
angle between (ab), (ac): 

(l,0\n-Jn ac |l, ) = ^!^ = -cosC- (62) 

slab -tt-ac 

Furthermore, the sine of the dihedral angle can also be directly extracted: 

(l,±l|n-Xc |1,0> =±-J=e^< sinC, (1, 0| |1,±1) = T -L e =F fa & sin^ c . (63) 



Let us now prove the formula (|33[) . Similarly to (I62|) . one has: 

E ba ■ Ad(g bc )E ca = (1, 0| n£ g bc n ca |1, 0), (64) 

by definition of the different group elements here involved. Then, write g bc in terms of the variables of twisted 
geometries, (|2"6")l : 

E ba ■ Ad{g bc )E ca = (1, 0| (n^ 1 n bc ) e e^ T * (n^ 1 n ca ) |1, 0). (65) 

The matrix elements of e = e~ nTy in the representation of spin j are: (j,m\e\j,m') — (— iy~ rn 8 m .- m '- Hence, 
introducing a resolution of the identity: 

E ba ■ Ad(g bc )E ca = M n ba n bc \l,m) (-^-"e*"*" (l,-m| n^ 1 n ca |1,0). (66) 

m=-l,0,l 

So (|33|) comes from explicitly writing down the three terms in the sum and the matrix elements. 

Quantization of the Hamiltonian in the large spin limit 

To evaluate the action of H c b l c on coherent spin network states, let us first collect a few results. The action of the 
generators of su(2) is what we expect from states with a nice semi-classical behaviour: 

f\j,n(N)) = (-i)jN\j,n(N))+o(j). (67) 
It is natural to see the 3- vector N in the Lie algebra su(2). Then, its components on the spherical basis (T TO ) m= _i j o,i 



are: 



N m = (l,m|n(7V)|l,0), m = -l,0,+l. (68) 
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SU(2) invariant operators acting on a coherent intcrtwiner naturally commute with the group action in (|4"5]) . In 
particular, these results can be used to compute the action of E ab ■ E ac . Each triad operator inserts a generator, so 
that we have to consider the action of T l ab ® T ac ) on the coherent intertwiner. It commutes with the group action, 
and finally leads to: 

(kTXc) s^"' n ^ ~ -j abJac (N ab ■ N ac ) S {J".».^}, (69) 

jabjac COSC S ^ (70) 

Acting with E ba ■ Ad(gb c )E ca is a bit more involved, but there are no conceptual difficulties. First the adjoint action 
of gbc is a Wigner matrix D^ 1 ' in the spin 1, which recouples to the matrix elements of g bc in the state, using 
Clebsch-Gordan coefficients: 

AB \9bc) ^kpVdbc) ~ 2^ °^ fcA+fc °B pB+p ^ A+k B+pV9bc) ■ [<!-) 
77=-l,0,l 

The question is then whether (and since the answer is yes, how) the coherent state \jb c , n bc ) in the state also receives a 
shift by r] = — 1,0, 1 consistently. This actually comes from the triad operators. Combining (|67l) and (|68|) . the action 
of E ba is recast so that the quantity of interest is a state in the tensor product of the representation spaces rlj ic 

\jbc,nbc) ® Tlba 1 1,0) = ^ ((1,H n bc n ba |1,0)J |jbc, «6c) ® «6c|l, m). (72) 
m=-l,0,l 

On the right hand side, we have just introduced the identity on Hi as: Ylm n bc\l> m )(l, m \ n b c • The quantity into 
brackets can then be evaluated thanks to the formula (|32|) which has already been used for twisted geometries. Quite 
clearly, this introduces the cosine of the dihedral angle 4> h ac when m = 0, and its sine, with some phase exp(=Fia£ a ) 
when m = ±1. Then, a key technical point which leads to the final expression is due to a careful inspection of the 
scaling properties of the Clebsch-gordan coefficients in Wj (g) Hi for large j: 

\j,n(N))(g>n(N)\l,m) ~ |j + m, n{N)). (73) 

This is exact when m = 1, but only holds asymptotically for m = —1,0. It enables to reintroduce a coherent state, 
with a shift on the spin, on the right hand side of ([721 : 

\jbc,n bc ) ® Ji6o|l,0) ^ -cos<^ c libc^bc) + -^e - * *- sin^ c |j bc + l,n fec ) - "^ eM " Q sin 0ac libc - l,"6c) (74) 

Obvisouly, a similar result holds for (jb c , n cb e\, on the other end of the link (6c), upon contraction with E ca . Gathering 
these different pieces, we get: 

(E ba ■ Ad(^)£ M ) S {3 -" n « l} ~ jabjac COS ^ COS0g a ^ C - 5 sin ^ c SU1 ^ (e-^ca+^J^c+l + e i(^ +0 s ^-l 

(75) 
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